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SUMMARY 


Analysis of a flexible beam with displaced end-located inertial masses is 
presented* The resulting three-dimensional mode shape is shown to consist of two 
one-plane bending modes and one torsional mode. These three components of the mode 
shapes are shown to be linear combinations of trigonometric and hyperbolic sine and 
cosine functions* Boundary conditions are derived to obtain nonlinear algebraic 
equations through kinematic coupling of the general solutions of the three governing 
partial differential equations* A method of solution which takes these boundary 
conditions into account is also presented. A computer program has been written to 
obtain unique solutions to the resulting nonlinear algebraic equations. This 
program, which calculates natural frequencies and three-dimensional mode shapes for 
any number of modes, is presented and discussed. 


INTRODUCTION 

With the advent of the Space Shuttle, a new class of spacecraft has been 
emerging in recent years. This class of spacecraft, which includes large antennas, 
platforms, space stations, etc., will require large lightweight space structures. 
These systems will be highly flexible, and will have a large number of significant 
vibrational modes which (unlike conventional rigid spacecraft) can no longer be 
ignored while designing control systems. Because of small inherent damping and a 
large number of elastic modes, these large flexible systems will need complex 
control laws in order to accomplish the required degree of precision in their 
pointing and maneuvering. These control laws, however, cannot be designed unless 
accurate vibrational characteristics of the system can be obtained. 

Many flexible space structures, as well as other terrestrial systems, can be 
represented by a single one-span beam with masses at both ends. Two shuttle based 
experiments which can be represented by this system are the Solar Array experiment, 
which flew in September 1984, and the Spacecraft Control Laboratory Experiment 
(SCOLE) (see ref. 1), which is a laboratory experiment in the design stage. For 
those structures whose end-mass centers of mass coincide with the principal axis of 
the beam, and have no product of inertia, simple one-plane vibrational analysis for 
each of the two planes (i.e., x-z and y-z planes) and one independent torsional 
vibration analysis are sufficient. (See ref. 2 for analysis of these types of 
systems.) However, if the centers of mass of the end masses do not lie on the beam 
axis, or if the products of inertia are not zero, then the three modes mentioned are 
coupled into one complex three-dimensional mode. 

This paper presents one method of approaching the problem of three-dimensional 
vibration analysis of a uniform beam with offset inertial masses at each end. The 
method assumes that the system is adequately described by three governing partial 
differential equations (i.e., two one-plane bending and one torsional bending 
partial differential equations). With proper boundary conditions, taking into 
account the center of mass displacements, and products of inertia, the natural 
frequency and mode shape can be obtained for any number of modes. This 
paper presents such a set of boundary conditions as well as a computer 
program that makes the solution readily obtainable for any given set of parameters. 
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SYMBOLS 


A cross sectional area of beam 

ApBi.CpDi coefficients of x-z plane mode shape equation 
A2,B2>C2,D2 coefficients of y-z plane mode shape equation 
A^,B3 coefficients of torsional mode shape equation 


c.m. 

center of mass 


El 

bending stiffness 

for beam when (EI)„ = (EI)„ 

x y 

<BI) X 

x-z plane bending 

stiffness 

(EDy 

y-z plane bending 

stiffness 


G modulus of rigidity 

Ip polar moment of inertia 

I xx ,Iyy,I zz x-, y-, and z-axis moments of inertia, respectively 

^xxO’^yyO’^zzO x ~» a n d z-axis moments of inertia, respectively, at 

z = 0 on beam 

^xxL^yyL’^zzL x “» and z ~ ax ^ s moments of inertia, respectively, at 

z = L on beam 


X xy 

xy 

product 

of 

inertia 



*xyO 

xy 

product 

of 

inertia at 

z = 0 

on beam 

X xyL 

xy 

product 

of 

inertia at 

z = L 

on beam 


L length of beam 


M x ,M y ,M z moments about x-, y-, and z-axes, respectively 


M x Q,M y Q,M z Q moments about x-, y-, and z-axes, 
z = 0 on beam 


respectively, at 


M x L,My L ,M z L moments about x-, y-, and z-axes, respectively, at 
z = L on beam 


mQ mass at z = 0 on beam 

m^ mass at z = L on beam 

p(t) common time solution of partial D.E's 

p x (t),py(t),p z (t) separate time solutions of xz-plane, yz-plane, and z-axis 

torsional P.D.E f s , respectively 
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r x (z),r x (e) xz-plane mode shape 
r y (z),ry(e) yz-plane mode shape 
t time 

u(z,t) beam displacement in xz-plane 

v(z,t) beam displacement in yz-plane 

v x> Vy shear forces in x- and y-directions, respectively 

V x o>VyQ shear forces in x- and y-directions, respectively, 

at z = 0 on beam 

V x L,VyL shear forces in x- and y-directions , respectively, 
at z = L on beam 

X,Y,Z position variables 

Z(uO eigenvalue matrix 

a phase angle (rad) 

B x ,8l mode shape variable for xz-plane 

6 y> 3 2 mode shape variable for yz-plane 

mode shape variable for z-axis torsion 

Axq , A yg c.m. displacements in x- and y-directions, respectively, 
at z = 0 on beam 

Ax^jAy^ c.m. displacements in x- and y-directions, respectively, 
at z = L on beam 

e dimensionless position variable (e = z/L) 

0(z), 0(e) z-axis torsional mode shape 
& K ,Q x (z y t) angular displacement about x-axis 
0y,0y(z,t) angular displacement about y-axis 
p density of beam 

<f>(z,t) angular displacement about z-axis 

aj natural frequency common to all three governing partial D.E ! s 

0 ) , 0 ) ,oj natural frequency of vibration of xz-plane, yz-plane, and z-axis 
x y z torsional bending modes, respectively 
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THE GOVERNING DIFFERENTIAL EQUATIONS 


The governing partial differential equations for the beam shown in figure 1 are 
comprised of two one-plane bending equations and one axial torsion equation. These 
differential equations all assume small displacements and slopes, uniform 
distribution of stiffness El and density p, and the torsional equation is derived 
specifically for a circular shaft. These three equations are the following (refs. 2 

and 3) : 

3 2 u(z,t) _ (EI) x 3 4 u(z,t) , 1N 

» A 9 z 4 0> 

for x-z plane bending, and 

.2 , (El) .4 . 

3 v(z , t) _ v y 3 v(z,t) 

3z 4 (> 

for y-z plane bending, and 


3%(z ,t) _ G 3 2 ^(z,t) 
3t 2 p 3z 2 


(3) 


for z-axis torsional bending* In equations (1), (2), (3), z denotes the 
independent space variable along the z-axis, t is the time, u, v, and <j> 
respectively denote the x- and y-axis bending and torsional displacement. (EI) x , 
(EI)^, G, P and A respectively denote the x-z and y-z plane bending 
stiffness, the modulus of rigidity, the density, and the cross-sectional area. 

These three equations are solved by separation of variables with the following 
substitutions (refs. 2 and 3): 


u(z,t) = r x (z)p x (t) 


(4) 


and 


v(z,t) = r y (z)p y (t) 


•1 

(5) 


and 


(j)(z,t) = e(z)p z (t) 


( 6 ) 
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Where r x , and 0 denote the bending and torsional mode shapes, and p x , Py, 

p 7 denote the corresponding functions of time* Substituting equations (4), (5), 
and (6) into equations (1), (2), and (3) and rearranging terms, equations (1), (2), 
and (3) are transformed, respectively, into: 


d 2 p x (t > 


dt 


P x (t) 


(EI) ? 

pA 


r (z) 
x 


d4r x (t) 

dz 4 


(7) 


d 2 Px (t) 

dt 2 


P x (t > 


(EI)^ 

pA 


r (z) 

y 


d 4 r y (z) 

dz 4 


( 8 ) 


and 


d 2 Pz (t) 


p z (t) 


dt 


G JL 

p 9(z) 


d 2 0(z) 

dz 2 


(9) 


Equations (7), (8), and (9) can be true if, and only if, both sides of each equation 
are equal to a constant* If the respective constants are chosen to be “w , -w , 
and -a) , the following six ordinary homogeneous differential equations are 
obtained* 

For the x-z plane bending: 


d p x (t> 


dt 


— + Vx (t) = 0 


( 10 ) 


and 


d4 r x (t) 4 

— 5— - 6 x r x (2) ■ ° 
dz 


do 


where 


4 _ pA 2 
e x _ ( EI ) X “ x 


( 12 ) 
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For the y-z plane bending: 


d p * (t) . 2 


dt 


2 + Vy (t) = 0 


and 


where 


d4r y ( z > 4 

2L g r (z) ■- 0 

dz 4 y y 


A _ oA 2 
e y " (EI) y w y 


( 13 ) 


(14) 


(15) 


For the z-axis torsional bending: 


d 2p z (t) 

2 — + V>z (t) = 0 

dt 2 2 


(16) 


and 


d 2 e(z) 

dz 2 


- 0 ), f 0(Z) = 0 
Z Cj 


(17) 


The system which is being modeled will be considered to be vibrating with the same 
frequency w in all three independent modes. This simplification gives all three 
modes the same time dependent governing equation: 

+ h) 2 p(t) = 0 . (18) 

dt 

[i.e., P x (t) = p y (t) = p z (t) = p(t)]. 
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The solution to equation (18) is (refs. 2 and 3): 


p(t) = cos (cot + a) (19) 


where a is a phase angle. 

The solutions to the position dependent equations (11), (14), and (17) with 

co = to = to = to are found to be (refs. 2 and 3): 

• x y z 

r x (z) = Aj sin B x z + cos B x z + Cj sinh B x z + Dj cosh g x z (20) 


r y (z) = A 2 sin gyZ + B£ cos ByZ + C 2 sinh B y z + D 2 cosh g y z (21) 


0 (z) = A 3 sin B z z + B 3 cos 8 z z ( 22 ) 


where 



2x1/4 
(0 J 



2a/4 
to J 



(23) 

(24) 

(25) 


Equations (20), (21), and (22) are more convenient to use when the position variable 
is transformed into a nondimensional form. For this reason the variable e = z/L, 
where L is the length of the beam, is used. After substitution, equations (20), 
( 21 ), and ( 22 ) become: 


r x ( e ) = Aj sin B^e + Bj cos Bje + sinh B^e + Dj cosh Bie 


(26) 


r y (e) = sin 02 e + ®2 cos ®2 e + C 2 s i n ^ &2 e + °2 cos ^ &2 e 
9(e) = A^ sin 63 ^ + B 3 cos B 3 G 


(27) 

(28) 
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where 


6 


e 


6 


4 

1 


4 

2 


2 

3 


pA 2 t 4 
w l 


pA 2 4 

tei r w L 

y 


£ 2 t 2 

a) L 


(29) 

(30) 

(31) 


Boundary Conditions for Three-Dimensional Vibrations of a Beam 
With Displaced Inertial End Masses 

The configuration being considered is a beam with inertial masses at both ends 
with x- and y-axis offsets* (See fig* 1*) The offset center of mass, along with 
the product of inertia, cause kinematic coupling between the x-z and y-z plane 
bending modes and the z-axis torsional mode* Figures 2 and 3 show the moment and 
shear force reactions being considered in the configuration. 

The following relationships between shear, moment, and beam displacement are 
used in the boundary conditions (ref. 4): 


V = -(El) 

X X 9z 3 


(32) 


,r _ / ft n 3 V(z,t) 
‘“’r 3 z 3 




- -(ei) 


y 3z 2 


y X 3z Z 


(33) 


(34) 


(35) 


M = GI M Z -— 
z p 3z 


(36) 
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When e = z/L is substituted into equations (32) through (36), the following 
relationships are obtained: 


(EI) X 

9 3 u(e,t) 

L 3 

8s 3 

1 

w 

M 

9 3 v(e,t) 

L 3 

3e 3 

1 

w 

M 

9 2 v(e,t) 

L 2 

9 e 2 

<EI) * 3 2 

u(e,t) 

L 2 

9 e 2 


9»(e,t) 
2 L 9e 


(37) 

(38) 

(39) 
(AO) 
(Al) 


Shear Forces at z = 0 

Referring to figure 2, the first boundary conditions involve the shears V x q 
and VyQ which are described by the following relationships: 


'xO 


- ( E I ) at z = 0 

X 9z J 


(A2) 


V. 


yO 


.(ED a . yfa.t) at z - 0 
y 9z j 


(A3) 


Setting equations (A2) and (A3) equal to the mass iiiq times the corresponding 
components of acceleration yields the following relationships. 




3 u(z , t) _ 3 <f>(z, t) 


at 


o 


at 


. (E I) 


(AA) 


3z^ 


ITK 


3 2 v(z,t) ^ 3 2 (J)(z,t) 

2 0 2 


3t 


3 1 


“ '(ED 


3 _v(z , t) 


x ^3 
3z 


(45) 
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The following substitutions will be used to transform equations (44) and (45) 
usable form: 


u( e, t) = r x (e)p(t) 


v(e,t) = r y (e)p(t) 


4 >(e,t) = 0 (e)p(t) 


where 


= z/L 


Using these four relationships, equations (44) and (45) become: 


2 

m o d " p( ~ 2 ~ [ r x (e) “ A yo 0(e U 

dt 


(EI) d 3 r (e) 

~ o p(t) 

L d e 


? at e = 0 . 


m 0 d ■ P( 2 -- [ r y < e) + A 0 9(e U = " 

dt 


(EI) y d J r y (s) 


de' 


p(t) 


J 


Equation (18) can be rewritten as: 


d 2 p(t) 
dt 2 _ 2 

W = 


Using equations (52), equations (50) and (51) can be rewritten as: 


2 (EI) x d r x (e) 

-w m o [ r x ( e ~ Ay 0 9 (e)] 3 3 — 

L de 


-03 m Q [r (e) + Ax Q 0(e)] = - 


(EDy Q r y (e) 


de' 


into 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 
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or after rearranging terms: 


d r x (e) m 0 4 

T“ = ^AL 6 J r x (e) + Ay O 0(e) ] at e = 0 


de~ 


( 55 ) 


where gj is given by equation (29). 


d r (e) m , 

1 = g 2 [r v (e) + AXq 0 (e)] at e 


def 


pAL M 2L y' 


= 0 


(56) 


where 62 is given by equation (30). 


Shear Forces at z - L 


Again referring to figure 2, the next set of boundary conditions involves the 
shears V xL and at z = L on the beam. These are derived with relationships 

similar to those used for the first two (V x q and V^q), ttie on ly difference being 
in the sign convention used in the derivation of the governing differential equa- 
tions (ref. 3)- For the end z = L the following relationships apply: 


and 


El, + \ 

__ 2 3 u(e,t) T 

V T “ n n 3t Z lj 

xL t 3 . 3 

L 9e 


El, »3 , * 

V = - at z = L 

yL L 3 9e 3 


(57) 


(58) 


Setting equations (57) and (58) equal to the mass times the corresponding 

component of acceleration, the following relationships are obtained: 


'xL “ 


'yL = \ 


9 u(e,t) _ 9 (J)(e,t) 

? a ^T 9 

dt L dt 


9 2 v(e,t) 3 2 j>(e,t) 

2 ax L .2 


(El) ^3 , > 

x 9 u(e,t) 

T 3 3 

L 9e 


/ at e = 1 


dt 


dt 


(EI) V 9 3 v(e,t) 

L 3 9e 3 ^ 


(59) 


(60) 
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Using relationships (46) through (49) and (52) in equations (59) and (60) and 
rearranging, the following boundary condition equations are obtained: 


d 3 r (e) til 

3 " ' = pAL M~ r x (e) + A y L 9(e) ] at e = 1 (61) 


where g^ is given by equation (29). 


d3r y (e) "Y 4 

~3 = ^AL e 2[" r y (e) " Ax L 0(e) ] at e = 1 
de 


where g 2 is given by equation (30). 


(62) 


Bending Moments at z = 0 

Referring to figure 3 the next two boundary conditions involve the moments M x q 

and M q at z = 0 on the beam# Equations (34) and (35) are equated with the 

following relationships (ref. 5) which ignore all nonlinear coupling and require 

that all products of inertia except I are zero: 

xy 


M = I 0+1 0 

x xx x xy y 


(63) 


M = I 0 + I 0 

y yy y xy x 


(64) 


Combining equations (63) and (64) with equations (34) and (35), respectively, and 
applying them to the end at z = 0 one obtains: 


M x0 = 1 n ® + 1 n® = “( EI ) -- v( § , - t) 

xu xxO x xyO y y 2 

9 z 


M 0 = I n e + I ft e = (EI) 
yu yyO y xyO x x 2 

dZ 


(65) 


at z = 0 


( 66 ) 
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The angular displacements 0 X and 0y can be approximated by the slopes of the y- 
and x-displacements, respectively: 


0 x (*,t) = “ 


9v(z, t) 
9z 


(67) 


e (z,t) 


9u(z,t) 

9z 


( 68 ) 


(The sign-convention used for rotations is the one corresponding to the standard 
right-handed coordinate system.) Using equations (67) and (68), equations (65) and 
(66) become: 


_ r 9 3 v(z,t) + j 9 3 u(z,t) = _/ EI \ 9 2 v(z,t) 
xx° 2 xyO 2 y - 2 


9z 9t 


9z 9t 


9z 


(69) 


> at z = 0 


3 3 2 

T 9 u(z,t) _ 9 y(z,t) _ , . 9 u(z,t) 

X yyO 2 xyO 2 ^ i; x . 2 

9z 9t 9z 9t 9z 


(70) 


J 


Using relationships (46), (47), and (49), equations (69) and (70) become: 


d 2 p(t) 

dt 2 


I n dr (e) I n dr (e) 
xxO y + xyO x 


de 


de 


(El) d"r (e) 

— 2 ~^ 

L de 


d 2 p(t) 

dt 2 


I n dr (e) 
yyO x 

L de 


I n dr (e) 
xyO y 

L de 


(EI) x d^r x (e) 

2 2 
L de 


p(t) 


(71) 


(72) 


Using equation (52) in equations (71) and (72), the following two boundary condition 
equations are obtained: 


d 2 r y (e) $2 


-I 


dr^Xe) dr x (e) 


de 


pAL 


xxO de 


+ I 


xyO de 


at e = 0 


(73) 
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where 02 given by equation (30). 


d 2 r x (e) 


de 


pAL 


dr (e) 

-I — + T 

yyO de 


dr y (e) 
xyO de 


at e - 0 


(74) 


where 0 ^ is given by equation (29). 

Bending Moments at z = L 

The fourth pair of boundary conditions involves the moments at the z = L end 
of the beam (see fig. 3). This set uses equations very similar to those used in 
set III. The only difference is in the sign convention necessary to satisfy the 
governing partial differential equations (ref. 3 ). The x- and y-moments at z * L 
are given by: 


M xL - (El) 


9 v(z,t) 
y dz 2 


at z = L 


MyL * -(EI)/J%Ei at t - L 
9z 


(75) 


(76) 


Setting equations (75) and (76) equal to equations (63) and (64), respectively, and 
using relationships (67) and ( 68 ), the following equations are obtained: 


3 3 v(z,t) T 3 3 u(z,t) _ a 3 2 v(z,t) - t 

^177?"" ‘“V"? - “ 


(77) 


3 3 2 

x 3 u(z,t) T 3 v(z , t ) _ /citN 3 u(z,t) _ _ _ T 

- - <EI) *-^“ at 2 - L 


(78) 


Using the relationships (46), (47), and (49), equations (77) and (78) become: 


d 2 p(t) 


dt 


I T dr (e) I T dr (e) 
xxL y xyL x 

L de L de 


(El) d r (e) 

__L 2— — p(t) at e = 1 (79) 

IT de 
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p(t) at e = 1 (80) 


d 




dr x (e) 

de 


I xyL 

L 


dr y ( e ) 

de 



d 2 r 


x 

de 


(e) 

2 


Using equation (52) in equations (79) and (80) and reducing gives rise to the 
following boundary condition equations: 


d 2 r y (e) 82 


de pAL 


dr (e) dr (e) 

I — X - I — ^ 

xxL de xyL de 


at e = 1 


where 82 is given by equation (30) • 


d 2 r x (e) bJ 


d e pAL [_ 


dr x (e) 


dr (e) 


- I 


yyL de xyL de 


at e = 1 


where 8 ^ is given by equation (29) « 


(81) 


(82) 


Torsional Moment at z = 0 

The ninth boundary condition involves the z-moment at z = 0 (see fig* 3). 

This moment is caused by the mass mQ and moment of inertia I 2z q according to the 
following relationship (refs. 3 and 5): 


m = 1 i_-fe.il + 

z0 L zz0 2 
dt 


m. 


& _ 9 2 u( z .t) Ay 
2 “0 , 2 fly 0 


dt 


dt 


(83) 


This moment is countered by the internal beam moment given by equation (36), 
Setting equation (83) equal to equation (36) gives the following: 


9 <|)(z,t) 

zzO . 2 
dt 


m. 


9 v(z,t) 


Ax. 


dt 


9 u( z ) t ) 
dt 2 



GI 


94>(z»t) 
p 9z 


(84) 


at z = 0 
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Using the substitutions given by equations (46) through (49), the following is 
obtained: 


d 2 p(t) r T 

2 l zzO 
dt 


0(e) + nip 


[r y ( e )Ax Q 


r x (e)Ay 0 ]} 


^*P 99(e) 
L 3 c 


p(t) 


( 85 ) 


at e = 0 

and rearranging terras, the following 

+ V y o r x (e) ] at e m 0 (86) 

where 0^ is given by equation (31). 


Using the substitution given by equation (52) 
boundary condition equation is formed: 


d9(e) 

de 


pLI 


't _I zzO 0(e) " m 0 Ax 0 r y (e) 


Torsion Moment at z = L 

The tenth and final boundary condition involves the moment at z = L (see 
fig. 3). This moment follows the same relationship given in equation (83) but the 
countering moment changes sign as in the other cases at z = L. This countering 
moment is given by: 


M T = -GI at z = L 

zL p 3z 


(87) 


The following equation is obtained, which is similar to equation (84); 


“L 


a 2 v(z,t) a 2 u(z,t) ... 

5 Ax = — A y T 

at 2 L at 2 L 


+ J 3 j>(z,t) = _ GI 3<f>(z ,t) 
zzL 2 p 3z 

o t 


( 88 ) 


at z = L 

Using the substitutions given by equations (46) through (49), equation (88) is 
changed to: 


d 2 p(t) , 
dt 2 1 ZZL 


GI 

9(e) + mjAx L r y ( e ) - Ay^e)]} = - 


99(e) 

9e 


p(t) 


(89) 
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Using equation (52) in equation (89) and rearranging terms gives the following: 


^dT 1 “ TErtW’fe) + ‘ Wx (c) I at E ■ 1 


( 90 ) 


where 63 is given by equation (31). 

By substituting equations (26), (27), and (28) and appropriate values of e 
(i.e. , e = 0 at boundary 1 and e = 1 at boundary 2) into equations (55), (56), 
(61), (62), (73), (74), (81), (82), ( 86 ), and (90), the following ten linear equa- 
tions are obtained, respectively: 




er 0 \ 


(91) 


_ ( 6 2 m 0 ) | 2 m 0 

A 2 \ P AL j a 2 + L 2 \ P AL /"2 ' \ pAL /" 3 


6 2 m 0 Ax 0l 
D„ + l — ~~ tt~ ~ B_ = 0 


(92) 


"Vl 


®l m L 

p A L Sln (6 1 ) ” COS^) 

A 1 + 

pAL cos(B]|) + s in (Bp 


B l m L 

pAL -smh(B 1 ) - cosh(B 1 ) 

C 1 + 

01% 

p - AL - - cosh(B 1 ) + sinh^) 

r 1 r 




p AL Sin(0 3 ) 


A 3 + 


N 

pAL cos(0 3 ) 


b 3 = 0 


(93) 




B 2 m L 

pAL sln(B 2 } - cos(B 2 ) 

A 2 + 

[pAL COs(B 2 ) + sln( e2 ) J 


B 2 m L 

— A ~ L ~ ■ sinh(B 2 ) - cosh($ 2 ) 


& 2 ,n T 

~pAL~ cosh ^ 2 ) + slnh(g 2 ) 


6 9 m T Ax t 


A 3 + 


Bn™. A x t 


b 3 - 0 


(94) 
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%4 A 1 * V + (%% - *2 + p¥ = 2 + °2 ' ° 




S 2 0 l I xvL 


|^—| {[cos(B 1 )]A 1 - [sin(B 1 )]B 1 + [coshC&p]^ + [sinh(B 1 )]D 1 } 



3, 


r 3-r 1 

+ 

bPxxl 

- cos(g 9 ) - sin( g_) 

- P AL 

A 2 + 

^2^xxL 

sin(g 2 ) - cos(g ) 



[ pAL 


$2 ^xxL ^2 ^xxL 

+ *— cosh(g ? ) - sinh(g 9 ) C„ + - *— sinh(g 9 ) - cosh(g 9 ) = 0 (97) 

pAI/ Z ' Z pAL J Z Z 


— - y ~ cos ( g. ) - sin( g. ) A + — -- 3X - sin( g. ) - cos(g. ) B. 

pAL pAL 


+ ^ cosh(g.) + sinh(g-) C. + - - -- - yy ~ sinh(g.) + cosh(g.) D. 

pAIT pAL 


e?M. 


3 ^— j { [cos(g 2 )]A 2 - [sin(g 2 )]B 2 + [cosh(g 2 )]C 2 + [sinh(g 2 ) ]D 2 } = 0 (98) 


e 3 m 0 Ay 0\ [ 6 3 ra 0 Ax 0\ 


+ A 3 + KeH B 3 - ° 
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Bo Ay, nu\ 

~~ p ~ L ~ i ~~j { [ slr| ( ) ]A l + [ cos( B 1 ) ]Bj + [sinh(B l )]C 1 + [cosh(B l )]D 1 } 


+ 


f % Ax L m L | 
PLI 


{ [-sin(B 2 ) ]A 2 + [cos(8 2 )]B 2 + 


[-sinh( 8 2 )]C 2 + [-cosh(B 2 ) ]D 2 } 


®3^zzL 

pLI 


sin(8 3 ) + cos(8 3 > 


A 3 + 


^3*zzL 

pLI 


cos( 8 3 ) - sin(6 3 ) 


b 3 = 0 


( 100 ) 


In equations (91) through (100), 8 1 , B 2 > and 63 are given by equations (29), 

(30), and (31), respectively. 


Obtaining Nontrivial Solutions 

Equations (91) through (100) can be written in vector-matrix form as follows: 


[Z(oj)] 


= 0 


( 101 ) 


Where Z(w) is the 10 x 10 coefficient matrix whose entries are functions of w 
(see equations (29), (30), and (31) for 8 j , 82 > and 63 )* Nonzero solutions 
(Ap Bp Cp Dp A 2> Bp Cp Dp Ap B^) exist only when the determinant of 
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Z(m) is zero. Therefore, the first step in obtaining nontrivial solutions is to 
obtain the real solutions of the nonlinear equations: 


det[Z( u )] = 0 (102) 

where det[ ] denotes the determinant. 

A solution 03 * is substituted back into equation (101) and a degenerate system 
(usually of rank 9) of algebraic equations is obtained. Choosing one coefficient 
and equating it to an arbitrary value (usually unity) the remaining nine 
coefficients can be uniquely determined for each solution w*. A computer program 
(BEAM3D), which calculates the nontrivial solutions, is discussed below. 


THE COMPUTER PROGRAM 

The computer program BEAM3D was written to obtain nontrivial solutions of 
equations (91) through (100) by the method discussed in the preceding section (see 
appendix A for a listing). BEAM3D was written assuming a symmetric cross section. 
Therefore, there is only one bending stiffness El = (EI) X * (EI)y* This was done 
simply because the equations are more accurate for symmetric or, more specifically, 
circular cross sections. This comes from the assumption used in deriving the 
governing partial differential equation for torsional vibration (eq. (3)) described 
in the first section. The only difference between the boundary condition equations 
used in the computer program and equations (91) through (100) is that 0^ = 02 
the program (see eqs. (29) and (30) for relationships of 0^ and B 2 ) since 
(EI) X = (EI) y . 

The boundary condition equations (91) through (100) are contained in ten 
separate subroutines named XSHRl, YSHR1, XSHR2, YSHR2, XM0M1, YM0M1 , XM0M2, YM0M2, 
ZM0M1 , and ZM0M2, respectively. The large number of trigonometric functions in 
these equations necessitated the use of additional variable names (for sin(g), 
cosh(g), etc.). 

The solution of equation (102) is computed in subroutine EIGEN. This sub- 
routine obtains values of the determinant of the 10 x 10 coefficient matrix cal- 
culated in the program external function FUNC. EIGEN checks for sign changes in the 
value of the determinant as well as changes in sign of the slope to find regions of 
possible roots. Once a region is found that contains a root, the root-solving sub- 
program SECBI is used to calculate the exact root. 

Subroutine SOLVE substitutes the root calculated in "EIGEN" into nine of the 
original ten boundary condition equations with one of the ten coefficients set equal 
to one to form a degenerate (rank 9) system. This system is then solved by using 
the system subprogram GELIM (for Gauss (- Seidel) elimination). 

Subroutine PHNORM divides the three mode shape equations by the appropriate 
factor (i.e., the square-root of the sura of the integrals of the mode shapes 
squared, integrated from zero to L over the space variable) in order to normalize 
them. This is useful for obtaining dynamic response to external forces or moments, 
or in control system studies. The only property of the mode shapes that changes 
because of normalization is the magnitude. 
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PROGRAM OPERATION 


Input* - A total of 22 input variables is needed to use BEAM3D. Use of a 
consistent set of units is required (meters-kilograms-seconds or feet-pounds- 
seconds). A description of each input variable can be found in the listing (see 
appendix A) in subroutine INNPUT. The only limitation on input is that at least one 
of the products of inertia ^ X y0 or *xyL must be nonzero to avoid a lower than 
rank 9 system. This is not really a limitation since the product of inertia 
couples the three bending modes in the first place. 

Output. - The output computed by BEAM3D includes the natural frequency of 
vibration, the normalized mode shape equations for x- and y- bending and z-axis 
torsion, and the plots corresponding to these mode shapes. A sample case giving 
both input and output is given in appendix B. 


CONCLUDING REMARKS 

A method of obtaining the natural vibration frequencies and mode shapes in 
three dimensions for a system comprised of a uniform beam with off-centered inertial 
masses at both extremities has been presented. The equations of motion were derived 
for this configuration taking into account the kinematic coupling resulting from the 
product of inertia and the offset end masses. The boundary conditions resulted 
in a set of nonlinear algebraic equations, the solutions of which yield the 
modal frequencies and mode shapes for any number of modes. A computer program 
was presented, which computes the modal frequencies and mode shapes for any 
desired number of modes. Since the mode shapes are comprised of trigonometric 
and hyperbolic sine and cosine functions, they can be readily differentiated 
to obtain the mode-slopes, which are required in control system studies. 
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APPENDIX A - COMPUTER PROGRAM BEAM3D 



o o o o o o o oooo 


I » 


c 

c 


10 


I J 
Iti 


PROGRAM BEAM3D ( INPUT* OUTPUT , TAPE? ) 

PROGRAM BEAM3D CALCULATES THE X-Z PLANE, Y-Z 
PLANE BENDING AND THE Z AXIS TORSION FOR A BEAM 
WITH INERTIAL MASSES WITH 'X* AND "V OFFSETS 
AT BOTH ENDS. 

COHMON/BEAM/ROAL, L,EI ,M1 ,M2, 1X1 , 1X8, IV1 , IY2, IXV 1 , IXV8 
1 , IZ1 , IZ8, ROLI , DX1 , DX8, DV1 , DV8, RHO, G, NMODE 
COMMON/TRIGOS/OMEGA1, OMEGA, BTA1,BTA2,SN,CS,SNH,CSH,SN8, CSS 
COMMON/SOLS/ A<1O),DET,DETSAUE,AMTRX1O<10,10>,AMTRX9<9,9>, 
1C0LUMN(9, 1 ),DETDIFF 

REAL 1X1, 1X2, IV1, IY2, IXV1,IXV2,IZ1, IZ2,L,M1,M2 
EXTERNAL FUNC 

PSEUDO INITIATES THE PLOTTING ROUTINE 
CALL PSEUDO 

SUBROUTINE INNPUT IS UHERE ALL INPUT UARIABLES 
ARE SIMPLY URITTEN IN AS THE TEXT 
CALL INNPUT 

OMEGA1 SETS THE STARTING POINT FOR THE NATURAL 
FREQUENCY OMEGA1 (RAD/SEC) 

OMEGA1 — .01 

DETSAUE AND DETDIFF ARE USED AS UORK UARIABLES 
IN SUBROUTINE EIGEN 
DETSAUE-1. 

DETDIFF-0. _ _ 

NATURAL FREQ'S. AND MODESHAPES ARE FOUND FOR THE 
FIRST NMODE NUMBER OF MODES 
DO 10 I -1, NMODE 
CALL EIGEN 

THE NEXT LINE PASSES UP THE TRIUIAL ZERO FREQUENCIES 
IF(BTA1.EQ.0. ) QO TO S 
CALL SOLUE 
CALL NORM (I ) 

CALL MPLOT 
CALL OUTPT(I) 

CONTINUE 

CALL C ALPLT ( 0 • , 0 • , 999 ) 

STOP 

END 
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SUBROUTINE EIGEN CALCULATES THE EIGENUALUES OF THE 
COEFFICIENT MATRIX 

SUBROUTINE EIGEN 

COMMON/TRI G0S/0MEGA1, ONEGA, BTA1 , BTA2, SN, CS, SNH, CSH, SN2, CS2 
COMMON/SOLS/ A ( 10 ) , DET, DETSAUE, AMTRX10 ( 10,10), AMTRX9 (9,9), 
1C0LUMN(9, 1 ),DETDIFF 

C0MM0N/BEAM/R0AL,L,EI,M1,M2,IX1, 1X2, IV1, IV2,IXV1, IXV2 
1, IZ1, IZ2,R0LI,DX1,DX2,DV1,DY2,RH0,G,NMCDE 
DIMENSION EPS(3 ) 

FUNC IS THE EXTERNAL THAT CALCULATES THE DETERMINANT 
OF THE 10 BY 10 COEFFICIENT MATRIX AMTRX10 
EXTERNAL FUNC 

REAL 1X1, 1X2, IY1, IY2, IXY1, IXY2,IZ1, IZ2,L,M1,M2 
EPS(l)-10**<-8.) 

EPS(2)-10**(-12) 

EPSO)-EPS(l) 

THIS IS UHERE THE FREQUENCY GUESSES ARE INCREMENTED 
5 0MEGA1 ■OMEGA1+0 * 01 
DET-FUNC(OMEGAl) 

THE FOLLOUING TUO IF-STATEMENTS LOCATE THE AREAS THAT 
MAY CONTAIN A ROOT (UHERE THE UALUE OF THE DET. CHANGES 
OR UHERE THE SIGN OF THE SLOPE CHANGES) 
IF(DETXDETSAUE.LE.0.) GO TO 10 
IF((DETDIFF*(DET-DETSAUE)).LE.0.) GO TO 10 
DETDIFF-DET-DETSAUE 
DETSAUE-DET 
GO TO 5 

10 OMEGA2*OMEGAl-.03 
DETDIFF-DET-DETSAUE 
DETSAUE'DET 

SECBI IS A ROOT FINDING ROUTINE IN FTNMLIB THAT USES 
A COMBINED SECANT-INUERSE QUADRATIC INTERPOLATION 
SAFEGUARDED BY BISECTION . (NOTE* THIS ROUTINE DOES 
NOT UORK UELL ON DOUBLE ROOTS, IF YOU HAUE DOUBLE 
ROOTS IT IS ADUISEABLE TO USE A DIFF. ROOT-SOLU. ROUTINE 
CALL SECBI ( OMEG A2 , OMEGA 1 , . 001 , FUNC , EPS , ROOT , IERR ) 



XF(XERR»EQ«3) URITE<7#10O) 

IF(XERR*EQ»8 ) URITE<?,200) 

IF(IERR*EQ«9) URITEC?#300) 

IF( IERR*EQ.8 ) GO TO 5 
IFC IERR.EQ.9 ) GO TO 5 

100 FORNAT( IX# 'SECBI TOOK MORE THAN 50 ITERATIONS TO FIND FREQ 1 ) 

200 F0RNAT(1X#'SECBI UAS GIUEN INPROPER INITIAL COND'S* SEE U0L.2 
1 C2*5 B ) 

300 FORNAT< IX# 'A ROOT UAS NOT FOUND IN THE INTERUAL GIUEN TO SECBI') 
ONEGA-ROOT 
CALL TRIG(ONEGA) 

RETURN 

END 
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SUBROUTINE PHNORM CALCULATES THE FACTOR NEEDED TO 
NORMALIZE THE MODE SHAPES 

SUBROUTINE PHN0RM(BETA*AC0*AC1*AC2*AC3*ELSQ*PHN*IM> 

B-BETA 

A1-AC1/ACO f A3-AC2/AC0 S A3-AC3/AC0 
S-SIN(B) S C-COS(B) S SH-SINHCB) S CH-COSHCB) 

S2 a SIN(2*XB ) S C2«C0S(2.XB) S SH2-SINHC2.XB > • CH2-C0SH(2.XB ) 
PHN* (B-S8/2. )/2.+AlXX2X(B+S2/2. )/2.+A2XX2X(SH2/2.-B )/2. 
1+A3XX2X ( SH2/2 . +B )/2 . +A1X < 1 . -C2 >/2 . +A2X ( SXCH-CXSH ) 
2+A3X(SXSH-CXCH+i. >+AlXA2X(CXCH+SXSH-l . )+AlXA3X(CXSH+SXCH ) 
3+A2XA3X ( CH2-1 « )/2. 

PHN-PHN/B 

PHN*SORT(PHN>XELSQXAC0 
URITE(?,201>PHN 
C COMPUTE APPROX. NORM USING SH-CH 

X«BX( 1 .+(Ai+A3 )XCAl-A2)+A3XX8)/2. 

! -S2X( 1 . +A1 >X< 1 .-Al >/4.+AiXC 1 .-C2 )/2.+A3-AlXA2-A2XA3/2. 
A23-A2+A3 

Y»SHXA23X(SX(A1+1. )+CX(Al-l. )+A23XSH/2. ) 

PHHl-X+Y 

PHN1 - SORT (PHN1/B )XELSQXACO 
WRITE (7*208 )PHN1 
IF(IM.GE.6)PHN-PHN1 

201 F0RMAT(1X*XN0RM USING EXACT F0RMULA-X*E12.5 > 

202 FORMAT ( IX* XMMAPPX NORM PHN1- X*E18.S) 

RETURN 

END 

C 

C 
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SUBROUTINE NORM DIVIDES ALL COEFFICIENTS BV THE 
FACTOR "PHN* WHICH IT RECEIVED FROM PHNORM 

SUBROUTINE NORN(N) 

COMMON/TRIGOS/ OMEGAl, OMEGA, BTA1,BTA3,SN,CS,SNH,CSH,SN2,CS2 
COMMON/SOLS/ A ( 10 > , DET, DETSAVE, AMTRXiO ( 10,10), AMTRX9 (9,9), 
1C0LUMNC9, 1 ),DETDIFF 

COMMON/BEAM/ROAL, L , El , Ml , MS, 1X1, 1X2, IV1, IV2, IXY1, IXY2 
1, IZ1, IZ2,R0LI,DX1,DX2,DY1,DY3,RH0,G,NM0DE 
REAL L 

ELSQ*SQRT(L) 

CALL PHNORMCBTAl, A( 1 ), A(S ), AO ),A<4 ),ELSG,XPHN,N ) 

CALL PHN0RM(BTA1, A(5 ),A(6 ),A(? ), A(8),ELSQ,YPHN,N ) 

B-A(10)/A(9) 

ZPHNSQ-A(9)**2*L/BTA3*((B**2+l)*BTA2/3.+CB**2-l)*SIN(3.*BTA2)/4 
l-B*C0S<3.*BTA2)/S,+B/2. ) 

PHN-SQRT ( XPHN**3+YPHN**3+ZPHNSQ ) 

DO 10 K-1,10 
A(K)-A(K)/PHN 
10 CONTINUE 
RETURN 
END 
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SUBROUTINE SOLUE CREATES AND SOLOES THE 9 BV 9 
REDUCED COEFFICIENT MATRIX OBTAINED BV CHOOSING 
AN ARBITRARY VARIABLE AND DROPPING ONE OF THE 
EQUATIONS 

SUBROUTINE SOLUE 

C0MM0N/BEAM/R0AL,L,EI,M1,M2, 1X1, 1X2, IY1, IY2, IXY1, IXY2 
1,IZ1,IZ2,R0LI,DX1,DX2,DY1,DY2,RH0,G,NM0DE 
COMMON/ SOLS/ A(10), DET, DETSAVE, AMTRX 10(10,10), AMTRX9 (9,9), 
1C0LUMN(9,1),DETDIFF 
DIMENSION IPIU0T(9),UK(9) 

REAL IX1,IX2,IV1,IV2, IXY1, IXV2, IZ1, IZ2,L,M1,M2 
NMAX,N,NRHS AND IFAC ARE ALL NEEDED IN GELIM 
NMAX-9 
N-9 

NRHS-1 
IFAC-0 
CALL S.ET(A) 

IARB DESIGNATES THE ARBITRARY VARIABLE 
I ARB "l 

ARB DESIGNATES THE VALUE GIVEN TO THE ARBITRARY VARIABLE 
ARB-1 

A( IARB ) --ARB 

ONE OF THE FOLLOUING 10 EQUATIONS IS COMMENTED 
SO AS TO BE IGNORED 
CALL XSHR1 ( COLUMN (1,1)) 

CALL XSHR2 ( COLUMN ( 1, 1 ) ) 

CALL VSHRK COLUMN (2,1 ) ) 

CALL YSHR2 ( COLUMN (3,1)) 

CALL XM0H1 ( COLUMN (4,1)) 

CALL XM0M2( COLUMN (5, 1 ) ) 

CALL VMOMK COLUMN (6, 1 ) ) 

CALL YM0M2 ( COLUMN (7,1)) 

CALL ZMONK COLUMN (8,1)) 

CALL ZM0M2( COLUMN (9, 1 ) ) 

C THE FOLLOUING DO-LOOP CALCULATES THE 9BY9 MATRIX 

21 DO 10 J-1,10 
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CALL SET(A) 

IFC J«EQ«IARB) GO TO 10 
K-J 

IFCJ.GT.IARB) IC-K-1 
AC J )-l • 

ONE OF THE FOLLOUING TEN EQUATIONS IS COMMENTED 
SO AS TO BE IGNORED 
CALL XSHR1CAMTRX9C1,K)> 

CALL XSHRS ( AMTRX9 C 1 , K > ) 

CALL YSHR1 CAMTRX9C3,IC > ) 

CALL YSHR8 ( AHTRX9 C 3, K ) ) 

CALL XMOM1 CAMTRX9C4,K > ) 

CALL XM0M2 ( AMTRX9 C 5, 1C ) ) 

CALL YMOM1 CAMTRX9C6,K > ) 

CALL YM0M2 ( AMTRX9 C ? , K ) ) 

CALL ZMOM1 C AMTRX9C8,K) ) 

CALL ZM0M2C AMTRX9C9,K > ) 

10 CONTINUE 

SUBROUTINE GELIM IS IN FTNMLIB AND SOLUES N BY N MATRICES 
CALL GELIMCNMAX,N, AMTRX9,NRHS, COLUMN, IPIVOT, IFAC,UK, IERR ) 
URITEC7,101) IERR 
101 FORMAT C3X, "IERR IS*,2X,I3) 

THE FOLLOUING LOOP ASSIGNS THE VALUES FOUND FOR THE 
9 NON-ARBITRARY VARIABLES FOUND IN GELIM IN ADDITION 
TO THE ONE ARBITRARY VARIABLE TO THE ORIGINAL ARRAY "AC 10) 


** 


DO 20 1-1,10 
IFC I »EQ* IARB ) GO TO 30 
IFC I *GT* IARB ) GO TO 40 
AC I )-COLUMNC I, 1 ) 

GO TO 20 
30 ACD-ARB 
GO TO 20 

40 AC I ) -COLUMN C 1-1, 1 ) 

20 CONTINUE 
RETURN 
END 


to 

VT3 
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SUBROUTINE MPLOT FORMS THE PLOTS OF THE THREE 
INDEPENDANT MODE SHAPES FOR EACH MODE 

SUBROUTINE MPLOT 

COMMON/SOLS/ A(10),DET, DETSAUE, AMTRX10C 10, 10 ) , AMTRX9 (9,9) 
1,C0LUMN(9, 1 ) 

COMMON/TRIGOS/OMEGA1 , OMEGA, BTA1, BTA2, SN, CS, SNH, CSH, SN2, CS2 
DIMENSION XLABC3 >, YLAB(3),XUER(3 ), YUERC3 ),T0RLAB(3 ),T0RUER<3 ) 
XLAB( 1 )«10HTHE XZ-PLA 
XLAB(2 )"10HNE MODE SH 
XLAB(3)»10HAPE (E) 

XUER( 1 )«10HX DISPLACE 
XUERC2)-10HMENT RX(E) 

XUER(3)»10H 
YLAB(1)«10HTHE Y2-PLA 
YLAB (2 )-10HNE MODE SH 
YLAB(3)-10HAPE (E) 

YUER( 1 )"10HY DISPLACE 
YUER(2)"10HMENT RY(E) 

YVER(3)»10H 

TORLAB ( 1 ) - 10HTORSIONAL 
TORLAB ( 2 ) • 10HMODE SHAPE 
TORLAB(3)-10H (E) 

TORUER ( 1 ) ■ 10HANGULAR DI 
TORUER ( 2 ) - 1 0HSPLACEMENT 
TORUER(3)*10H (E) 

CALL CRUNCH (A(l ),A(2),A(3),A(4),XLAB,XUER,BTA1 ) 

CALL CRUNCH ( ACS ),A(6),A(7),A(8), YLAB, YUER,BTA1 ) 

CALL CRUNCH (AC9),A<10),0.,0., TORLAB , TORUER, BTA2 ) 

RETURN 

END 
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SUBROUTINE CRUNCH USES INFOPLT TO GENERATE PLOTS 
FOR MPLOT OF EACH INDEPENDANT SUB-MODE 

SUBROUTINE CRUNCH < A, B, C, D, HLABEL, VLABEL, BTA ) 

DIMENSION HLABEL C 3 ) , ULABEL (3 )*BTA(1 ), EPSLN ( 1010 ), R C 1010 ) 

JJ-1 

DLTX IS THE INCREMENT ADDED TO EPSLN TO OBTAIN 
POINTS FOR THE PLOT OF THE MODE SHAPE 
DLTX- *001 

IEC,N,KX,KY,XMIN,XMAX,YMIN,YNAX,PCTPTS,NXMC,NYMC, 

ISYM,SX,SY,XOFF AND YOFF ARE ALL VARIABLES NEEDED 
IN INFOPLT* 

IEC-1 

N-1001 

KX-1 

KY-1 

XMIN-0.0 
XMAX-1 . 

YMIN— *03 

YMAX-.03 

PCTPTS-0.00 

NXMC-30 

NYMC-30 

ISYM-0 

SX-7. 

SY-5* 

X0FF-.75 
YOFF -.75 

EPSLN(JJ) IS THE HORIZONTAL COMPONENT OF THE PLOT 
UHERE AS R(K ) IS THE VERTICAL* 

EPSLNC 1 )-0. 

DO 30 K-1#N 

R(K)-A*SIN(BTA(1 )*EPSLN< JJ) >+B*COS(BTA< 1 )*EPSLN( JJ ) >+C*SINH(BTA< 1 ) 
1 *EPSLN ( J J > >+D*COSH ( BTA ( 1 )*EPSLN C J J ) ) 



UEPSIN ( J J ) )+D*COSH ( BTA ( i >*EPSIN ( J J ) ) 

5 EPSLN( JJ+1 )«EPSIN( JJ )+DLTX 

jj.jj+1 
20 CONTINUE 

CALt INFOPLT< IEC # N # EP$tN( 1 )«KX»R( 1 )#KY#XI1IN # XHAX # YNIN,VI1AX# 
lPCTPTS,NXnC,HLABEL,NYMC,ULABEl,ISVM,SX,SV,XOFF,YOFF> 

RETURN 

END 
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SUBROUTINE TRIG FINDS VALUES OF THE TRIGONOMETRIC 
FUNCTIONS. IT DOES THIS ONCE FOR EACH VALUE OF ONEGA. 

SUBROUTINE TRIG(DLTA) 

COMMON/BEAM/ROAL,L,EI,Mi,M2,IXl,IXa,IYl,IY2,IXYi,IXY2 
1, IZ1, IZ2 # R0LI,DX1,DX2,DY1,DY2*RH0,G > NM0DE 
COMNON/TRIGOS/OMEGA1 , ONEGA# BTA1 # BTA2# SN* CS# SNH # CSH* SN2, CS2 
REAL IXl,IX2,IYi,IY2,IXYl,IXY2,IZl,IZ2,L,Mi,MB 
BTA1 • ( R0AL/EI*DLTA**2*L**3 )** . 25 
BTA2»DLTA*L$ ( RHO/G )tt • 5 
SN-SIN(BTAl) 

CS-COS(BTAl) 

SNH-SINH(BTAl) 

CSH-COSH(BTAl) 

SN2-SIN(BTA2) 

CS2-COS(BTA2) 

RETURN 

END 




c*i 
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FUNC CALCULATES THE DETERMINANT OF THE 10 BY 10 MATRIX 
FUNCTION FUNC (BUMS) 

C0MM0N/BEAM/R0AL,L,EI,M1,M2,IX1, IX8,IY1, IY2,IXY1, IXY2 
1 , IZ1, IZ2, ROLI , DX1 , DX3, DY1 , DY2, RHO, G, NMODE 
COMMON/SOLS/ A<10>,DET,DETSAUE,AHTRX10(10,10>,AMTRX9<9,9> 
1, COLUMN (9,1 ) 

DIMENSION B ( 10 ) , IPIUOT ( 10 ) , UK (10 ) 

REAL IX1,IX2,IY1,IY8,IXV1,IXY2,IZ1,IZ8,L,M1,M8 
NMAX-10 
N-10 
IDET-2 

CALL TRIG (DUNS ) 

DO 10 I«l,10 
CALL SET ( A ) 

ACD-l. 

CALL XSHR1(B(1)> 

CALL XSHR2(B(2) ) 

CALL YSHR1(B(3) ) 

CALL YSHR8(B(4) ) 

CALL XM0M1(B(5>> 

CALL XM0M8(B(6) ) 

CALL YH0M1(B(?>) 

CALL YM0M8(B(8) ) 

CALL ZM0N1(B(9) ) 

CALL ZH0M8(B(10)) 

DO 80 J-1,10 
AMTRX10( J,I )*B( J) 

20 CONTINUE 
10 CONTINUE 

C DETFAC IS A DETERMINANT FINDING ROUTINE 

CALL DETFAC(NMAX,N,AHTRX10,IPIU0T,IDET,DUM,ISCALE,UK,IERR) 
FUNC ■ DUN* ( 10** ( 100BISCALE ) ) 

C URITE(?,282) DUM2,FUNC,ISCALE 
288 FORMAT (2X,E14«?,4X,E14*?,4X, 12) 

RETURN 

END 
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SUBROUTINE SET IS USED TO SET ALL THE COEFFICIENTS 
OF THE NODE SHAPE EQUATIONS EQUAL TO ZERO 

SUBROUTINE SETCA) 

DIMENSION A( 10) 

REAL IX1,IX8,IY1,IY8,IXY1,IXYB,IZ1,IZ8,L,M1,N8 

DO 10 1-1,10 

A(I>-0. 

10 CONTINUE 
RETURN 
END 
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SUBROUTINE XSHRl CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE SHEAR FORCE IN THE X-DIRECTION AT 
THE END Z-9. 

SUBROUTINE XSHRl (UX1> 

C0MN0N/BEAM/R0AL,L,EI,H1,M3, 1X1, 1X2, IY1, IY2, IXY1, IXY8 
1, IZ1,IZ2,R0LI,DX1,DX3,DY1,DY8,RH0,G,NM0DE 
CONMON/TRIGOS/ONEGA1 , ONEGA, BTA1 , BTA3, SN, CS, SNH, CSH, SN2, C$2 
COMMON/SOLS/ A ( 10 ) , DET, DETSAUE, AMTRX10 ( 10, 10 ) , AMTRX9 (9,9), 

1 COLUMN (9,1 ),DETDIFF 

REAL IXl,IX2,IYl,IY8,IXYl,IXY8,IZl,IZ2,L,Nl,Ma 
UX1— BTA1**4*(-A<10)*DY1+A<4)+A<2>)*N1/ROAL+AC3)*BTA1**3-A<1>*BTA 
1 1**3 

RETURN 
END 


SUBROUTINE XSHR8 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE SHEAR FORCE IN THE X-DIRECTION AT 
THE END Z-L. 

SUBROUTINE XSHRS(UXS) 

CONNON/SOLS/ A<lO>,DET,DETSAUE,ANTRX10<i0,lO),AHTRX9<9,9), 
1C0LUMN(9,1),DETDIFF 

C0NM0N/BEAM/R0AL,L,EI,N1,M2,IX1, IXS, IY1, IY2,IXY1,IXY2 
1 , IZ1 , IZ2, ROLI , DX1 , DXS, DY1 , DY8, RHO, G, NHODE 
C0NM0N/TRIG0S/0NEGA1 , ONEGA, BTA1 , BTA8, SN, CS, SNH, CSH, SN2, CS8 
REAL IX1,IX8,IY1,IY8,IXY1,IXV2,IZ1,IZ8,L,H1,N8 

VX2*( (A(4)XBTA1XX3XR0AL+A(3)XBTA1*X4XH2)XSNH-A(9)XBTA1XX4XDY2XM2 

1 *SN2+ < A ( 8 )*BTA1**3*R0AL+A ( 1 )*BTA1**4*N8 )*SN+ ( A ( 3 JXBTA1XX3XCSH- 

2 A( 1 )«BTA1*«3XCS )*ROAL+(-A< 19 )*BTA1**4*CS2*DY2+A(4 JXBTA1XX4XCSH+ 

3 A(2)*BTA1**4*CS)*N8)/R0AL 
RETURN 

END 
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SUBROUTINE VSHR1 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE SHEAR FORCE IN THE Y-DIRECTION AT 
THE END Z-e. 


SUBROUTINE YSHRl(UYl) 

CONNON/SOLS/ A<10),DET,DETSAUE,ANTRX10<iO,18),AHTRX9<9,9), 
1C0LUHN(9,1),DETDIFF 

C0HN0N/BEAH/R0AL,L,EI,Hl,H8,IXl,IX2,IVl,IYa,IXYl,IXY2 
1,IZ1, IZ2,R0LI,DX1,DX2,DY1,DYS,RH0,G,NN0DE 
CONNON/TRIGOS/ONEGAl , ONEGA, BTA1, BTAS, SN, CS, SNH, CSH, $N2, CS8 
REAL 1X1, 1X2, IY1,IY8, IXV1,IXV8, IZ1,IZ2,L,N1,H2 

UY1— BTA1**4*CA(10)*DX1+A(8)4A<6))*N1/ROAL+A(7)*BTA1**3-A<S)*BTA1 


1 **3 

RETURN 
END 


SUBROUTINE VSHR8 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE SHEAR FORCE IN THE Y-DIRECTION AT 
THE END Z-L. 


SUBROUTINE YSHR2(UY2> 

CONNON/SOLS/ A ( 10 ) , DET, DETSAUE, ANTRX10 < 10,10), ANTRX9 (9,9), 
1C0LUNN(9,1 ),DETDIFF 

C0HN0N/BEAN/R0AL,L,EI,H1,H8,IX1,IX8,IY1,IY2,IXY1,IXY2 
1, IZ1, IZ2,R0LI,DX1,DX8,DY1,DY2,RH0,G,NH0DE 
COHNON/TRIGOS/OREGA1 , ONEGA, BTA1 ,BTA8,SN,CS, SNH, CSH, SN8, CS2 
REAL IX1,IX2, IY1,IY2,IXY1,IXY8,IZ1,IZ8,L,N1,N8 

UY2- ( ( A ( 8 )*BTA1**3*R0AL+A ( 7 )*BTA1**4*H2 )*SNH+A < 9 )*BTAlX*4XDXa*H2 

1 XSN2+ ( A ( 6 )*BTA1**3*R0AL+A ( 5 )XBTA1**4*N2 )*SN+ ( A ( 7 )*BTA1X*3*CSH- 

2 A ( 5 )*BTA1**3*CS )*ROAL+ ( A ( 10 )*BTAi**4*CS8*DX2+A ( 8 )*BTA1**4*CSH+ 

3 A(6)*BTA1**4*CS)*N2)/R0AL 
RETURN 

END 


C>1 
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SUBROUTINE XN0N1 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE NONENT ABOUT THE X-AXIS AT Z-0. 

subroutine xnoni(nxi) 

connon/sols/ a ( 10 ) , det, detsaue, antrxio < 10, 10 ) , ahtrxs (9,9), 

1C0LUNNC9, 1 ),DETDXFF 

conhon/bean/roal,l,ei,ni,R2,ixi,ix2,iyi,iy2,ixyi,ixy2 
1, IZ1, IZS,R0LI,DXl,DX2,DYl,DY2,RH0,G,Nn0DE 
CONNON/TRIGOS/OHEGAl , ONEGA, BTA1 , BTA2, SN, CS, SNH, CSH, SN2, CS2 
REAL 1X1, 1X2, IY1, IY2,IXY1,IXY2, IZ1, IZ2,L,N1,N2,NX1 
NX1- BTA1**4*((A(7)*BTA1+A(5)*BTA1)*IX1-(A(3)*BTA1+A(1)*BTA1)*I 
1 XY1)/(L**2*R0AL)+A(8)*BTA1**2-A(6)*BTA1X*2 

RETURN 
END 


SUBROUTINE XN0N2 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE NONENT ABOUT THE X-AXIS AT Z-L. 

SUBROUTINE XN0N2CNX2) 

CONNON/SOIS/ A ( 10 ) , DET, DETSAUE, ANTRXIO ( 10,10), ANTRX9 (9,9), 
1C0LUNN(9,1 ),DETDIFF 

C0NN0N/BEAN/R0AL,L,EI,N1,N2, 1X1, 1X2, IY1, IY2,IXY1,IXY2 
1, IZ1, IZ2,R0LI,DX1,DX2,DY1,DY2,RH0,G,NN0DE 
CONNON/TR I GOS/ONEGA 1 , ONEGA , BTA1 , BTA2 , SN , CS , SNH , CSH , SN2 , CS2 
REAL 1X1, IX2,IY1, IY2, IXY1,IXY2, IZ1, IZ2,L,N1,H2,NX2 
NX2» BTA1X*4*(IXY2*(A(4)*BTA1*SNH-A(2)*BTA1*SN+A(3)*BTA1*CSH+A(1) 

1 *BTA1*CS)-IX2*(A(8)*BTA1*SNH-A(6)*BTA1*SN+A(7)*BTA1*CSH+A(5)* 

2 BTA1XCS) )/( LXX2XR0AL )+A ( 7 )*BTA1**2*SNH-A ( 5 )*BTA1**2*SN+A ( 8 XBTA 

3 1X*2*CSH-A(6)*BTA1**2*CS 
RETURN 

END 
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SUBROUTINE YM0M1 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE MOMENT ABOUT THE Y-AXIS AT 2-0. 


SUBROUTINE YMOMl(MYl) 

COMMON/SOLS/ A(10),DET,DETSAUE,AMTRX10(10,10),AMTRX9<9,9), 
1C0LUMN(9,1 ),DETDIFF 

C0MR0N/BEAM>'R0AL,L,EI,M1,M2,IX1,IX2, IY1,IY2, IXY1, IXY2 
1, IZ1, IZ2,R0LI,DX1,DX2,DY1,DY2,RH0,G,NM0DE 
C0MM0N/TRIG0S/0HEGA1 , OMEGA, BTAi , BTA2, SN, CS, SNH, CSH, SN2, CS2 
REAL 1X1, IX2,IY1, IY2, IXY1, IXY2,IZ1, IZ2,L,H1,M2,MY1 
MYl«*BTAi8848((A<3)8BTAl4-A<l)8BTAl)8IYi-<A<?)8BTAl+A(5)8BTAl)8IX 
1 Y1 )/( L8828R0AL )+A ( 4 )8BTA1882-A ( 2 )8BTA1882 

RETURN 
END 


SUBROUTINE YM0M2 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE MOMENT ABOUT THE Y-AXIS AT Z-l. 


SUBROUTINE YM0M2CMY2) 

COMMON/SOLS/ A ( 10 ) , DET, DETSAUE, ARTRX10 ( 10,10), AMTRX9 ( 9 , 9 ), 
1C0LUMN(9,1),DETDIFF 

C0MM0N/BEA!VR0AL,L,EI,R1,M2,IX1,IX2, IY1, IY2, IXY1, IXY2 
1, IZi, IZ2,R0LI,DXi,DX2,DYl,DY2,RH0,G,NM0DE 
COFiMON/TRIGOS/OREQAl , OMEGA, BTA1, BTA2, SN, CS, SNH, CSH, SN2, CS2 
REAL 1X1, 1X2, IY1, IY2, IXY1, IXY2, IZ1, IZ2,L,H1,H2,RY2 

RY2 0 B7A18848<IXY28(A<8)SBTA18SNH-A<6)88TAiSSN+A<7)8BTA18CSH+A(5) 

1 aBTAl«CS)-IY28(A(4)8BTAl*SNH-A(2)aBTAl«SN+A<3)»BTAl«CSH+A(l)*B 

2 TA1SCS) )/(L8828R0AL)+A(3)8BTA18828SNH-A(l )8BTA18828SN+A(4 )8BTA1 

3 8828CSH-AC2 )8BTA18828CS 
RETURN 

END 
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SUBROUTINE ZN0M1 CONTAINS THE BOUNDARY CONDITION 
DESCRIBING THE MOMENT ABOUT THE Z-AXIS AT Z-0. 

SUBROUTINE ZHOMKNZl) 

COnnON/SOLS/ A<lO>,DET,DETSAVE,AHTRXlO<lO,10),AriTRX9<9,9), 
1C0LUMN(9,1 >,DETDIFF 

coMMON/BEAM/ROAL,L,Ei,iu,M8,ixi,ix8,iYi,iYa,ixYi,ixYa 
l,IZl,IZ8,R0LI,DXl,DXa,DYl,DYa,RH0,G,NM0DE 
C0MM0N/TRIG0S/0NEGA1 , OMEGA, BTA1 , BTA2, SN, CS, SNH, CSH, SN3, CS2 
REAL IXl,IX8,IYl,IY3,IXYl,IXY2,IZl,IZa,L,ni,M8,MZl 
MZ1 - BTA8**8*<<AC8HA<6>>*DX1*M1/R0LI-(A<4>+A(S>)XDY1XM1/R0LI+A<1 
1 0>*IZ1/R0LI)+A(9>*BTA8 

RETURN 
END 


SUBROUTINE ZM0M3 CONTAINS THE BOUNDARY CONDITIONS 
DESCRIBING THE MOMENT ABOUT THE Z-AXIS AT Z-L. 

SUBROUTINE ZMOMS(MZS) 

CONNON/SOLS/ A<1O),DET,DETSAVE,ANTRX1O(10,10>,AMTRX9<9,9>, 

1 COLUMN (9,1), DETDIFF 

C0HM0N/BEAH/R0AL,L,EI,Hl,Ra,IXl,IXa,IYl,IY8,IXYl,IXYa 
1 , IZ1 , IZS, ROLI , DX1 , DX2, DY1 , DY2, RHO, G, NMODE 
COnnON/TRIGOS/ONEGAl , ONEGA, BTA1 , BTA8, SN, CS, SNH, CSH, SN8, CSS 
REAL IX1,IXS,IY1,IY3,IXY1,IXV3, IZ1,IZ8,L,N1,N8,NZ8 
NZ8— ( ( A ( 7 )*BTA8«8*DX8-A ( 3 )SBTA8X*8*DY8 XM8XSNH+ ( A ( 10 ) 3BTA8XR0LI+ 
1 A ( 9 >*BTA8**8*IZ2 XSN8+ ( A ( 5 XBTA8SX8XDX8-A ( 1 XBTASXXSXDYS XMSSSN-A 
8 ( 9 >*BTA8*CS8*R0LI+ ( ( A ( 8 >*BTA8**3XCSH+A ( 6 XBTA8XX8XCS )*DX8+ < -A ( 4 ) 

3 *BTA8**8*CSH-A(8 >XBTA8**3*CS )XDY8 >*M8+A< 10 )XBTA8*S8XCS8XIZ8 )/R 

4 OLI 
RETURN 
END 


i 
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SUBROUTINE INNPUT 

c0MM0N/BEAH/R0AL,L,Ei,Mi,H2,ixi,ixa,iYi,iYa,ixYi,ixYa 
1,IZ1,IZ8,R0LI,DX1,DX3,DY1,DY3,RH0,G,NM0DE 
REAL IXl,IXa,IYl,IY3,IXYl,IXYB,IZl,IZ2,L,Hl,Ha 
C "ROAL* IS THE NASS OF THE BEAN ALONE* 

ROAL-13.43 

0 , L" IS THE LENGTH OF THE BEAN. 

L-130. 

C "EX* IS THE BENDING STIFFNESS OF THE BEAN (SYMMETRICAL) 

EI-4. *10**7. 

C ■Ml" IS THE MASS OF THE LUMPED MASS AT Z-0. ON THE BEAM 

Ml -6366. 46 

C "MB* IS THE HASS OF THE LUMPED MASS AT Z-L ON THE BEAM 

M2-12.4a 

C 1X1 AND 1X3 ARE THE MOMENTS OF INERTIA ABOUT THE 

C X-AXIS AT Z-0. AND L RESPECTIVELY 

1X1 -905443.0 
1X8-18000. 

IY1 AND IY3 ARE THE MOMENTS OF INERTIA ABOUT YHE 
Y-AXIS AT Z-0. AND L RESPECTIVELY 
IV1-6789100. 

IY2*9336* 

IXYl AND IXY2 ARE THE PRODUCTS OF INERTIA AT Z-0,L RESP. 
IXY1-0.0 
IXY2*-?570f 

IZl AND iza ARE THE MOMENTS OF INERTIA ABOUT THE 
Z-AXIS AT Z-0.,L RESPECTIVELY 
IZ1-7086601. 

IZ2”27407# 

C 'PI* is THE POLAR MOMENT OF INERTIA ABOUT THE Z-AXIS 

PI*#1738 

C ‘“RHO 1 IS THE DENSITY OF THE BEAM. 

RHO-. 9089/PI 
ROLI-RHO*L*PI 

C DX1 AND DY1 ARE THE X AND Y DISPLACEMENTS OF THE POINT 



to 


C MASS AT 2*0* DX2 AND DV2 ARE THE SANE FOR Z»L. 

DX1-0.0 

DV1«0.0 

DX2-18.75 

DV2»-32*5 

C *G * IS THE NODULUS OF RXDQIDITV 

G-4.E+7/PI 

0 NNODE SPECIFIES THE NUNBER OF NODES TO BE SOLOED FOR 

NN0DE*5 

URITE(7,180) R0AL,L,EI,ni,N2,XXl, XX2,IY1, IY2,IXY1,IXY2, IZ1,IZ2 
1 , ROLI , DX1 , DX2, DY1 , DY2,RH0,G 
1C0 FORNATdX, 'THE BEAN CHARACTERISTICS ARE*,/, 

11X, *NASS»*,E14.7,/,1X, "LENGTH* 1 ,E14.7,/, IX, • STIFFNESS* ",E14. 7, /, 
21X, "I11*",E14.7,/,1X, "M2*",E14.?,/, IX, "IX1*",E14.7,/, IX, *IX2 a ",El 
34.7,/, IX, * IY1* " ,E14.7,/, IX, * IY2*" ,E14.7,/, IX, 1 IXY1* 1 ,E14.7,/,1X, 
4" IXY2* *,E14.7,/, IX, "XZ1*",E14.?,/,1X, "IZ2*" ,E14.7,/, IX, "RH8L8IP*" 

5, E14.?,/,lX,"DXi*",E14.7,/,lX,,"DX2*",E14.7,/,lX,"DYl*",E14.7,/,lX 

6, "DY8*",E14.7,/,lX,"DENSITY*",Ei4.7,/,lX,"G*",E14.?> 

RETURN 

END 
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SUBROUTINE OUTPT(I) 

C0MM0N/BEAH/R0AL,L,EI,Ml,Ha,IXl,IX2,IYl,IY2,IXYl,IXY2 
1, IZ1, IZ2,R0LI,DX1,DX2,DY1,DY2,RH0,G,NM0DE 
C0MM0N/TRIG0S/0MEGA1, OMEGA, BTA1,BTA3,SN,CS,SNH, OSH, SN2,CS2 
COMMON/SOLS/ Ad0),DET,DETSAUE,AMTRX10dO,10),AMTRX9<9,9> 

1. COLUMN (9,1) 

REAL IXi, 1X2, IV1, XV8, XXVI, IXY2, XZ1 , IZ2,L,M1,M2 
0MEGAHZ*0MEGA/(2.X3« 1415926 ) 

URITE(7,10) I,OMEGAHZ 

10 FORMATdX, a THE SOLUTION FOR MODE • , I2,///,2X, 'THE FREQUENCY OF 
1UIBRATI0N IS* ,E14.7 ) 

URITE(7,20) (A(J),J*1,4) 

20 FORMATdX, 'THE XZ-PLANE MODE SHAPE IS «,/,2X,E14.7, *XSIN(BETAUX/L 
l)+ , ,E14.7, , XC0S<BETAl*Z/L>+ - ,/,2X,E14.?, , *SINH(BETAl*Z/L>+ , ,E14*7, 
B'XCOSHdETAUZ/L)-) 

URITE(7,30 ) (A( J),J*5,8) 

30 FORMATdX, *THE YZ-PLANE MODE SHAPE IS • ,//,2X,E14*7, -XSINCBETAIXX/ 
1L)+ , ,E14.7, , *C0S(BETA1XZ/L)+",/, 

22X,E14.7, , *SINH<BETA1*Z/L>+\E14.7, , *C0SH<BETA1*Z/L) , > 

URITE(7,40) (A( J), J"9,10) 

40 FORMATdX, "THE TORSIONAL MODE SHAPE IS , ,//,2X,E14.7, "*SIN<BETA2*Z 
l/L)+ , ,E14♦7, , *C0S(BETA2*Z/L) , ) 

URITE(7,50) BTA1,BTA2 

50 F0RMATdX, , BETAl* , ,E14.7,/,lX, , BETA2* , ,E14.7) 

RETURN 

END 





APPENDIX B - SAMPLE CASE USING BEAM3D 


Figure B-l shows the Spacecraft Control Lab Experiment or SCOLE (ref. 1) 
configuration, which is the system to be analyzed in this test case. The following 
list contains all the parameters of the SCOLE geometry needed by the computer program 
BEAM3D: 


Mass of Space Shuttle = mg “ 6366.46 slugs 
Mass of Reflector = - 12.42 slugs 

Length of Beam = L = 130 feet 


Inertias of Shuttle 
at the Attachment Point 


Inertias of Reflector 
at the Attachment Point 


I I xx0 = 905,443 slug-ft 2 

I yyQ = 6,789,100 slug-ft 2 

2 

I zz0 = 7,086,601 slug-ft 
LyO = 0 slu S' ft2 

/ X xxL = 18 > 000 slu g~ ft2 
Z yyL = 9 » 336 slug-ft 2 

LzL = 27 » 407 slug-ft 2 
T X y L = 7,570 slug-ft 2 


Shuttle CM Location 


Reflector CM Location 


AXq = 0. ft. 

Ay 0 = °* ft ' 

Ax = 18.75 ft. 

L 

Ay = 32.5 ft. 

Li 


Material Properties 


pA = .09554 slugs/ft 

EI = 4. x 10 7 lb-ft 2 

pi = .9089 slug-ft 
P 

GI = 4. x 10 7 lb-ft 2 
P 


44 



SPACECRAFT CONTROL LAB 
EXPERIMENT (SCOLE) 




Figure B-l : Drawing of the SCOLE, Shuttle/ Antenna Configuration 

modelled for the Sample Case 



THE SOLUTION FOR NODE 1 


THE FREQUENCV OF VIBRATION IS *3?40493E+0O HZ* 

THE XZ-PLANE NODE SHAPE IS 

* 1616907E+00XSIN(BETA1XZ/L)+ -.1964167E+00*COS(BETA1*Z/L)+ 
-*1688419E+00*SINH(BETA1XZ/L)+ .1958759E+00*COSH<BETA1*Z/L) 

THE VZ-PLANE NODE SHAPE IS 

-*3957009E-01XSIN(BETA1XZ/ t)+ *6906827E-O1XCOS<BETA1*Z/L )+ 

»5848914E-01XSINH(BETA1XZ/L>+ -.6890796E-O1*COSH<BETA1*Z/L ) 

THE TORSIONAL NODE SHAPE IS 

-*3199889E-01XSIN(BETA2XZ/L)+ .158U62E-04*C0S(BETA2*Z/L> 

BETA1- .1192552E+01 
BETA2- .3374271E-01 


Figure B-2a: Natural Frequency and Mode Shapes calculated 

by BEAM3D for Mode #1. 



X DISPLACEMENT 



Figure B-2b : Projection of the first mode shape, calculated using BEAM3D, 

onto the xz-plane where the displacement r x is plotted 
versus the nondimensional position variable e. 




DISPLACEMENT 



ANGULAR DISPLACEMENT 0(e) 




TORSIONAL MODE SHAPE (e) 


Figure B-2d: The torsional deflection of the first mode shape, calculated 
— ' ‘ ‘ using BEAM3D, where 6 (in rad.) is plotted versus the 

nondimensional position variable £. 


THE SOLUTION FOR NODE 8 


THE FREQUENCY OF VIBRATION IS *3229025E+00 HZ* 

THE XZ-PLANE NODE SHAPE IS 

• 7265880E-01XSIN ( BETA1XZ/L )+ -.8441782E-01XCOSCBETA1XZ/D+ 
-*75062O1E-01XSINH<BETA1XZ/L>+ *8419519E-01XCOSHCBETA1XZ/L) 

THE YZ-PLANE NODE SHAPE IS 

*1857088E+O0SSIN(BETA1XZ/ L)+ -.1965677E+O0XCOS(BETA1XZ/L >+ 
- • 1676755E+00XSINH ( BETA1 XZ/L )+ . 1961255E+00XCOSHCBETA1XZ/L > 

THE TORSIONAL NODE SHAPE IS 

. 2S99605E-O2XSIN C BETA2XZ/L )+ - . 109020OE-OSXCOS ( BETA2XZ/L ) 

BETA1- .1294490E+01 
BETA2- .3975783E-01 


Figure B-3a ; Natural Frequency and Mode Shapes calculated 
by BEAM3D for Mode 92. 



X DISPLACEMENT 



Figure Projection of the second mode shape onto the xz-plane where the 

i displacement r x is plotted versus the nondime nsional position 

variable e. 


Y DISPLACEMENT 



ANGULAR DISPLACEMENT 6(e) 



THE SOLUTION FOR NODE 3 


THE FREQUENCY OF VIBRATION IS • 7487723E+00 HZ » 

THE XZ-PLANE NODE SHAPE IS 

•2297884E-01*SXN(BETA1XZ/L)+ -♦S911988E-01*COS<BETA1»Z/LH 
-»2345591E-018SINH(BETA1SZ/D+ ♦5907393E-O1*COSHCBETA1*Z/L) 

THE YZ-PLANE NODE SHAPE IS 

.2543413E-0i*SIN(BETAl*Z/ LH *.3262716E-02*COSCBETA1*Z/L>+ 
-*2523513E-01XSINH(BETAltZ/L)+ -.3312861E^02*COSH(BETAl*Z/L > 

THE TORSIONAL NODE SHAPE IS 

; 7256298E-01BSIN ( BETA2BZ/L H - ♦ 13123O7E-04BCOS ( BETA2BZ/L ) 

BETA1 ■ • 1971232E+01 

BETAS* .9219366E-01 


Figure B-4a ; Natural Frequency and Mode Shapes,: calculated 
by BEAM3D for Mode #3. 



X DISPLACEMENT r x (e) 


/ 



Figure B-4b: Projection of the third mode shape onto the xz-plane where the 

displacement r x is plotted versus the nondimensional position 
variable 


Y DISPLACEMENT 



ANGULAR DISPLACEMENT 6(e) 


Ln 






THE SOLUTION FOR NODE 4 


THE FREQUENCY OF VIBRATION IS *1244013E+01 HZ* 

THE XZ-PLANE NODE SHAPE IS 

• 6875555E-01XSIN ( BETA1SZ/L )+ -.638S680E-01*CO$(BETA1*Z/L>+ 

- « 6899687E-01 XSINH ( BETA1SZ/L )+ *63?5084E-01«COSH(BETAUZ/l) 

THE YZ-PLANE NODE SHAPE IS 

- • 1050247E+00SSIN ( BETA1XZ/ L)+ *940O61OE-01*COS<BETA1*Z/L)* 

* 1 O76795E+00XSINH ( BETA1XZ/L )+ - * 9384279E-01XCOSH ( BETAlXZ/l ) 

THE TORSIONAL NODE SHAPE IS 

. 1 139082E-01XSIN ( BETA2XZ/L >♦ - . 1239939E-05XCOS ( BETA2XZ/L ) 

BETA1- *2540828E+01 
BETAS- .1531709E+00 


Figure B-5a ; Natural Frequency and Mode Shapes calculated 
by BEAM3D for Mode // 4 . 
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Figure B-5b: Projection of the fourth mode shape onto the xz-plane where the 

displacement r x is plotted versus the nondimensional position 

variable e. 
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THE SOLUTION FOR NODE 5 


THE FREQUENCY OF UIBRATXON 1$ *2051804E+81 HZ* 

THE XZ-PLANE NODE SHAPE IS 

• 9723593E-01 tSIN ( BET AlftZ/L )+ -.9O63536E-01*COSCBETA1*Z/L)+ 
-•9?39?13E-01*SINH(BETA1*Z/L>+ .9O519O0E-01*COSHCBETA1*Z/L) 

THE VZ-PLANE NODE SHAPE IS 

*5767110E-8USIN<BETA1*Z/ L)+ -.5459663E-01*COS<BETA1*Z/L )+ 
-.583991?E-01*SINH(BETA1*Z/L)+ *5452724E-01*COSH<BETA1*Z/L) 

THE TORSIONAL NODE SHAPE IS 

- . 4658896E- 03XS I N ( BET A2BZ/L ) + • 3074803E-07SCOS ( BETA2XZ/L ) 

BETA1- «3263103E+01 
BETA2* « 25263 13E+00 


Figure B-6a : Natural Frequency and Mode Shapes calculated 

by BEAM3D for Mode #5. 
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Figure B-6b: Projection of the fifth mode shape onto the xz-plane where the 

displacement r x is plotted versus the nondime ns ional position 
variable e. 
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displacement Vy is plotted versus the nondimensional position 
variable e. 
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03 



Figure B-6d ; The torsional deflection of the fifth mode shape where 0 (in 

rad.) is plotted versus the nondime ns ional position variable e. 


o 
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